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1 Introduction 

In current theoretical physics there is a relevant number of theoretical investigations that 
lead us to believe that at the Big-Bang first moments, the geometry was not commutative 
and the dominating physics at that time was ruled by the laws of noncommutative (NC) 
geometry. Therefore, the idea is that the physics of the early moments can be constructed 
based on these concepts. 

The first published steps through this knowledge were given by Snyder [1] who believes 
that NC principles could make the quantum field theory infinities disappear. However, it 
was not accomplished [2] and Snyder's ideas were put to sleep for a long time. The main 
modern motivations that rekindle the investigation about NC field theories came from 
string theory and quantum gravity [3]. 

In the context of quantum mechanics for example, R. Banerjee [4] discussed how NC 
structures appear in planar quantum mechanics providing a useful way for obtaining them. 
The analysis was based on the NC algebra in planar quantum mechanics that was originated 
from 't Hooft's analysis on dissipation and quantization [5]. 
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It is opportune to mention here that this noncommutativity in the context of string 
theory mentioned above could be ehminated constructing a mechanical system which re- 
produces the string classical dynamics [6]. NC field theories have been studied intensively 
in many branches of physics [7]-[14]. 

In a very interesting paper, a parallel investigation was developed by Duval and 
Horvathy [15], where it was obtained the anomalous commutation relations for the co- 
ordinates obtained through the "Peierls substitution" [16]. From first principles, without 
using such unphysical limit, the authors introduced NC (quantum) mechanics starting with 
group theory and applied it to condensed matter physics, e.g., the Hall effect. The respec- 
tive Lagrangian approach was discussed in detail in subsequent papers [17]. Dunne, Jackiw 
and Trugenberger [18] justify the Peierls rule by considering the m ^ limit, reducing 
the classical phase space from four to two dimensions, parameterized by NC coordinates 
X and Y, whereas the potential becomes an effective Hamiltonian. 

In [19] the authors analyzed the IR and UV divergences and verified that Planck's 
constant enters via loop expansion. Here, differently, we make a non-perturbative approach 
and we will see that Planck's constant enters naturally in the theory via Moyal-Weyl 
product. 

A general algebra a-deformation of classical observables that introduces a general NC 
quantum mechanics was constructed in [20]. This a-deformation is equivalent to some 
general transformation for the usual quantum phase space variables. In other words, the 
authors discuss the passage from classical mechanics to quantum mechanics. Then to NC 
quantum mechanics, which allows to obtain the associated NC classical mechanics. This 
is possible since quantum mechanics is naturally interpreted as a NC (matrix) symplectic 
geometry [21]. 

In [22], the author constructed an extension of the well known Doplicher-Fredenhagen- 
Roberts NC algebra introducing the formalism which is now called in the literature as 
the Doplicher-Predenhagen-Roberts-Amorim algebra. In this formalism the NC parameter 
(6) is an ordinary coordinate of the spacetime and therefore it has a canonical conjugate 
momentum (vr). An extended Hilbert space was constructed together with all the ingredi- 
ents of a new NC quantum mechanics. But notice that both preserves the underlying NC 
relation [x^,x'^] = iO^^'^. For details, the interested reader can consult [23]. 

Back to our main subject here, in few words we can say that one way to obtain NC 
versions for field theories one have to replace the usual product of fields into the action by 
the Moyal-Weyl product, defined as 



where 6^'^ is a real and antisymmetric constant matrix. As a consequence, NC theories 
are highly nonlocal. We also note a basic NC property that the Moyal-Weyl product of 
two fields inside the action is the same as the usual product, considering that we discard 
boundary terms. Thus, the noncommutativity affects just the vertices. 

Some years back, [24] three of us have proposed a new formalism to generalize the 
quantization by deformation introduced in [20] in order to explore, with a new insight. 
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how the NC geometry can be introduced into a (commutative) field theory. To accomphsh 
this, a systematic way to introduce NC geometry into commutative systems, based on 
Faddeev-Jackiw symplectic formahsm and Moyal-Weyl product, was presented [25]. 

One important arena where NC ideas may play an important role is cosmology. If 
superstrings is the correct theory to unify all the interactions in nature, it must have played 
the dominant role at very early stages of our Universe. At that time, all the canonical 
variables and corresponding momenta describing our Universe should have followed a NC 
algebra. Inspired by these ideas some researchers have considered such NC models in 
quantum cosmology [26-28]. It is also possible that some residual NC contribution may 
have survived in later stages of our Universe. Based on these ideas some researchers have 
proposed some NC models in classical cosmology in order to explain some intriguing results 
observed by WMAP. Such as a running spectral index of the scalar fluctuations and an 
anomalously low quadrupole of CMB angular power spectrum. Among such proposals we 
may mention the following ones [29-33]. Another relevant application of the NC ideas 
in classical cosmology is the attempt to explain the present accelerated expansion of our 
Universe [34-36]. 

We have organized this paper as follows. In section 2, we introduce the generalized 
quantization by deformation formalism assuming a generic classical symplectic structure. 
We will construct a new star product which appears at first sight to suffer from the non- 
associative property disease. We will demonstrate that we can recover this crucial property 
for the star product. In section 3, we construct the NC versions of several Friedmann- 
Robertson- Walker cosmological models. These models may have positive, negative or zero 
spatial sections curvatures, they are coupled to different types of perfect fluids and they 
may have a positive, negative or zero cosmological constant. We have used the Schutz 
variational formalism in order to write the Hamiltonian of the models. In section 4, we 
perform a complete study of the evolution of the universes described by each model. We 
solve the dynamical equations for the scale factor. Based on our results, we conclude how 
the NC parameter modifies the evolution of the universe in different models. Also in this 
section, we present some special cases where our conclusions can be clearly verified. In 
section 5, based on our previous results, we write a function (A) that depends on three 
parameters that generalizes the cosmological constant. We discuss the possible scenarios for 
the Universe evolution predicted by A, depending on the value of a. The final considerations 
and the conclusions are depicted in the last section. 

2 The NC Generalized Symplectic Formalism 

The quantization by deformation [37] consists in the substitution of the canonical quanti- 
zation process by the algebra Ah of quantum observables generated by the same classical 
one obeying Moyal-Weyl product, i.e., the canonical quantization 





with ( = {qi,Pi), is replaced by the ?i-star deformation of Aq, given by 
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{h,g}r, = h-kn.g - g-khh , (2.2) 

where 

{h^gm = exp{^;^x^,,9f^^)4^)}/i(Ci)5(C2)la=C2=C > (2-3) 
with a,b = 1,2, ... , 2N and with the following classical symplectic structure 

c^a6=^_°/g'^ With i,j = 1,2,..., N , (2.4) 

which satisfies the relation 

^^'^W = 5c • (2.5) 

In the next section we will describe the conceptual basis that support our method and 
our cherished results. However, we think that it is crucial for the reader to understand 
mathematically the underlying equations used in this work. We will show that the crucial 
property of associativity is not lost in our formalism. 



2.1 Deformation quantization 

Notice that we have a subtle conceptual analogy between both star products, i.e., the 
one defined by the Moyal-Weyl product in Eq. (1.1) and the S-star deformation product 
defined in Eq. (2.3). In general the star products so defined are associative if the parameter 
is constant. If we do not have associativity, the product so defined is useless in physics. 
We will talk with more detail about this from now on. 

Kontsevich proved in [38] that any finite-dimensional Poisson manifold can be canon- 
ically quantized (deformation quantization). We will write now a quite brief introduction 
to deformation quantization following [38]. 

Let us define an algebra A = T{X,Ox) over TZ of smooth functions on a finite- 
dimensions C°°-manifold X. We construct a star product on A defined as being an asso- 
ciative i?[?i]-linear product on A[h]. This star product between / and g (f, g & A C A[h]) 
can be written as 



{f,9)^f*g = fg + hBi{f,g) + h^B2{f,g) +...eA[h] (2.6) 

where h \s a. constant parameter and Bi are bidifferential operators. A bidifferential op- 
erator can be understood as bilinear maps which are differential operator [38]. And the 
product of arbitrary elements of A[f}\ can be defined following the framework written in 
(2.6) and it is defined by the condition of linearity over TZ[fi\ 

( E * ( E = E ^ 91 n"^' + E ^"^'^^ ■ 

n>0 n>0 A:,Z>0 k,l>0,m>l 

For more details the interested reader can look in [38]. 
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The simplest of a deformation quantization is the Moyal-Weyl product for the Poisson 
structure on TZ"^ with constant coefficients, 

f*9 + fg + hY,^' dj{g) + Y E ^Ma) + ■ ■ ■ 

i,j i,j,k,l 
oo , jj n n 71 

E n«""(na.O(^>^(n».o<''> 

n=0 ji,... ,j„ fc=l /c=l fc=l 

where a*-' = —a^^. 

Let a = ■ ot^ di A dj he a Poisson bracket with variable coefficients in an open 
domain of TZ'^ [38]. Namely, a*-' is not a constant but a function of coordinates. Then the 
following star product gives an associative product modulo 0{h^) [38], 

f*9 = fg + hY,a'^diif)djig) + y j;i,j,fc,/a*^a^'a,afc(/)5,9z(<7) 
+ y E cx''dj{a'''){d.dk{f)di{g) - dk{f)dMg)) + 0{h') . (2.8) 

For us to demonstrate the associativity up to the second order means that for any three 
functions f,g and h we have that 

{f*g)*h = f*{g*h) (2.9) 

We can see clearly form (2.8) that the fact of a*-* being a non-constant parameter brings 
a different order for the Moyal-Weyl product. A simple comparison between Eqs. (2.7) 
and (2.8) can make us see that if we can not write a star product like the one described in 
(2.8) so we do not have associativity. We will turn back to this issue in a few moments. 

The quantization by deformation can be generalized assuming a generic classical sym- 
plectic structure S"*. In this way the internal law will be characterized by h and by another 
deformation parameter (or more). As a consequence, the S-star deformation of the algebra 
becomes 

{hMj^gm = exp{^;iS,feaf^^)4^)}/i(Ci)5(C2)lG=C2=C : (2-10) 

with a, 6= 1,2, ...,27V. 

This new star-product generalizes the algebra among the symplectic variables in the 
following way 

{h,g}hT: = ihYlah ■ (2.11) 

Notice that the new star product in (2.10) is defined in an analogous way as the Moyal- 
Weyl product. However, in Eq. (2.11) we see that the parameter is not constant. Hence, 
we can realize that from Eqs. (2.10) and (2.11) the associativity property was lost in Eq. 
(2.10). 

However, if we consider only terms up to h? and due to the smallness of h we can 
recover the associativity property of Eq. (2.10). So, the general expression given in Eq. 
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(2.10) hides the fact that the expansion of the exponential is physicahy valid only for 
terms proportional to 1, ^ and h"^ and therefore it is not completely correct. Our objective 
in Eq. (2.10) was to write the E-star product in a compact way. But we have to make this 
observation. 

On the other hand, we will see soon that Eq. (2.10) is not the starting point of our 
procedure. The main point is the generalized Dirac quantization defined in the next section. 

The motivation to construct a S-star deformed product like the one in Eq. (2.10) is 
to show that we can have in physics an alternative and at the same time associative star 
product with variable coefficients without losing the associativity property as we showed 
above, i.e., where the expansion stops in the second order of h. 

2.2 The NC approach 

In [20, 21], the authors proposed a quantization process to transform the NC classical 
mechanics into the NC quantum mechanics, through generalized Dirac quantization, 

{h,9}j: = ^^at^ \[Oh,0,h ■ (2.12) 

The relation above can also be obtained through a particular transformation onto the usual 
classical phase space, namely, 

C = TabC' , (2.13) 

where the transformation matrix is 

T=(l'' ~fA , (2.14) 

\ 12 Pv / 

with 0ij and /3jj being antisymmetric matrices. As a consequence, the original Hamiltonian 
becomes 

nca) no . (2.15) 

The corresponding symplectic structure is 

^ab=(_/l , (2.16) 

\ Ptj ) 

'^ij ~ ~T&\^ikPkj + Pik^^kj]- Due to this, the commutator relations look like 

[q'i^q']] = ^^ij ' 

[q'i^P]] =iK5ij + (yij) , (2.17) 

At this point, it is important to notice that a Lagrangian formulation was not given so 
far. Now, we propose a new systematic way to obtain a NC Lagrangian description for a 
commutative system. In order to achieve our objective, the symplectic structure S^fe must 
firstly be fixed and after that, the inverse of must be computed. As a consequence, an 
interesting problem arise: if there are some constants ( Casimir invariants) in the system. 
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the symplectic structure has a zero-mode, given by the gradient of these Casimir invariants. 
Hence, it is not possible to compute the inverse of However, in Ref. [39] this kind of 
problem was solved. On the other hand, if T,ab is nonsingular, its inverse can be obtained 
solving the next relation 



J^abix, y) J^'^y, z)dy = 5l5{x - z) , (2.18) 

which generates a set of differential equations since S'^^ is an unknown two-form symplectic 
tensor obtained from the following first-order Lagrangian 

C = A^,C-yiO , (2.19) 

as being 

5Ar (x) SA(-> (x) 

Due to this, the one-form symplectic tensor, Aq{x), can be computed and subsequently, 
the Lagrangian description, Eq. (2.19), is obtained also. 

In order to compute A^^(x), the Eq. (2.18) and Eq. (2.20) will be used, which generates 
the following set of differential equations 

9ijBjk{x, y) + {5ij + aij) Ajkix, y) = 6ik6{x - y) , 

Ajk{x, y)9ji + {6ij + aij) Cjk{x, y) = , 

- {6ij + aij) Bjk{x, y) + (3ij Ajkix, y) = , (2.21) 

Akjix,y){6ji + aji) + l3ijCjkix,y) = 6ik6{x - y) , 

where 



6Ag,{x) 


5Ag,^{x) 




Wj{y) 


5Ap,{x) 


M,,(x) 


Kiy) 


^v'jiy) 


5Ap>{x) 


5Ap,^ (x) 




5p'j{y) 



Bjk{x,y) = 
Ajkix, y) = 

C,k{x,y) = [ - ^-^^^ ) . (2.22) 

From the set of differential equations in Eq. (2.21), and the equations above, Eq. (2.22), 
we compute the quantities A(^i^{x). 

As a consequence, the first-order Lagrangian can be written as 

C = A^,Ca - ViO . (2.23) 

Notice that, despite (2.19) and (2.23) have the same form, in (2.23) the A(^'^ are completely 
computed through the solution of the system (2.21). In both we have a NC version of the 
theory as a consequence of the deformation in (2.14) and its corresponding symplectic 
structure in (2.16). 
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3 The noncommutative Friedmann-Robertson- Walker models 

Priedmann-Robertson- Walker (FRW) cosmological models are characterized by the scale 
factor a{t) and have the following line element, 

ds^ = -Nitfdt" + a{tf (y^^ + r^dn^^ , (3.1) 

where is the line element of the two-dimensional sphere with miitary radius, N{t) is 
the lapse function and k gives the type of constant curvature of the spatial sections. The 
curvature is positive for = 1, negative for k = —1 and zero for = 0. Here, we are 
using the natural unit system, where c = BvrG = 1. We assume that the matter content 
of the model is represented by a perfect fluid with four- velocity U^^ = 5q m. the co- moving 
coordinate system used, plus a cosmological constant (A) which can be either positive, 
negative or zero. The total energy-momentum tensor is given by. 



7)ti/ = {p + p)UyUy - pg^,y - Ag^y , (3.2) 

where p and p are the energy density and pressure of the fluid, respectively. Here, we 
assume that p = ap, which is the equation of state for a perfect fluid. 

Einstein's equations for the metric (3.1) and the energy momentum tensor (3.2) are 
equivalent to the Hamilton equations generated by the following super-Hamiltonian con- 
straint [40], 

n = ^-3ka + Aa^ + Pra"^" , (3.3) 

12a 

where Pa and Pt are the momenta canonically conjugated to a and T respectively, the 
latter being the canonical variable associated to the fluid. This super-Hamiltonian was 
derived using the Schutz variational formalism [41, 42]. The starting point to derive the 
NC version of the above cosmological models is the super-Hamiltonian constraint (3.3). 

In order to obtain a NC version for the FRW model, we apply the procedure described 
in the previous section. Initially, we have to write the zeroth-iterative Lagrangian of the 
system, which can be done directly from the super-Hamiltonian in (3.3), 

=Paa + /¥T-y(a,pa,r,PT), (3.4) 
where V = NVl is the symplectic potential and 

_p2 

= ^-3A:a + Aa^ + PTa~^". (3.5) 

12a ^ ' 

Notice that the system is treated classically via Symplectic Formalism. From now on we 
will follow the steps described in the last section. 

In the Lagrangian (3.4), the symplectic variables are identified easily as 



C = {a,Pa,T,PT,N) , 
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and the corresponding zeroth-iterative one-form canonical momenta 



AO) _ p AO) 



A 



(0) 







A 



(0) 








A. 



(0) 
T 



Pt 



(3.6) 



Calculating the one-form canonical momenta using the symplectic matrix definition, 

_dA^_dA^ 

we obtain directly the zeroth-iteration symplectic matrix, given by 

/O -1 0\ 
1 
0-10 
10 
VO 0/ 



(3.7) 



(3.8) 



However, this matrix is singular, which assumes the existence of constraints in the system, 
and it has the following zero-mode 



1 

Multiplying this zero- mode by the symplectic potential gradient we have that, 

4 



E 

1=1 



dV 



(3.9) 



(3.10) 



where is a constraint. This constraint can be introduced into the kinetic sector of the 
zeroth-iterative Lagrangian C^'^\ through the Lagrangian multiplier /?. In this way, the 
first-iterative Lagrangian can be written as 



= p^a + PtT + ni3-Nn 



(3.11) 



and the new set of symplectic variables is ^* = {a, Pa,T, Pt, N, 13). The corresponding 
first-iterative one-form canonical momenta A^i(S^^Y^'^ given by 



4'^ = Pa A 



A 



(1) 







A 



(1) 

Pa 
(1) 








(1) 

T 

(1) 



Pt 

n. 



(3.12) 



Thus, from the relation (3.7) we obtain the first-iterative symplectic matrix 



/ 



(1) 








-1 











dn 

da 


\ 




1 














dn 

dPa 






1 



























1 








dn 
dPr 


























v 


da 


On 

dPa 





dn 
dPr 








/ 



(3.13) 
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However, this matrix is also singular. If we multiply its zero-mode by the gradient of the 
symplectic potential we will find the same constraint obtained previously, 



E 



fJ'i 



dV 



n. 



(3.14) 



This result leads us to conclude that the system has a gauge symmetry. In accordance 
with the symplectic method, this symmetry must be introduced into the Lagrangian in 
Eq. (3.11) through the Lagrange multiplier S. So, the second-iterative Lagrangian can be 
written as 

£(2) = p^a + PtT + St) - A^O , (3.15) 

where = N — 1, i.e., the lapse function is equal to one, which is equivalent to the choice 
of a physical time, and the new set of symplectic variables is ^* = {a, Pa,T, Pt, N,r]). 
Using the symplectic matrix definition in Eq. (3.7) again, we have a non-singular second- 
iterative symplectic matrix. After a straightforward calculation, we obtain the inverse of 
the symplectic matrix, 

/ 1 



\ 
-1 
1 
-1 
-1 
V 1 / 



(3.16) 



It is important to remember that the elements of this matrix corresponds to the Poisson 
brackets among the symplectic variables, (/~^)*-' = {C The standout moment of the 
NC introduction through the method described in the last section lies in the assumption 
that the following relations for the brackets Poisson are. 



{a,T} = e 

{Pa,PT]=P. 



(3.17) 



These brackets are justified by the fact that T and its conjugated momentum Pt can be 
defined in terms of specific entropy and the potential e [40], which can be written as a 
function of energy and the position coordinate. Furthermore, in NC classical mechanics 
the brackets among the coordinates can be nontrivial. 

With all these values in mind, using the inverse of the symplectic matrix, including 
the NC brackets in Eq. (3.17), we can determine the symplectic matrix directly 



/ 





' 1 


-P 













-1 





-e 








1 


/3 





1 








139-1 


9 


-1 






















(/ 


39-1) 




. 








{1-P9) 






(3.18) 
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where ^50— 1 7^ 0, is a constraint. To proceed with the method, we use the symplectic matrix 
elements (3.18) and the relations in Eq. (3.7). The result is a set of partial differential 
equations, 



dAp^ 


dAa 


1 


da 


dPa P 


9-1 


BAt 


dAa _ 


-/3 


da 


dT (3. 


9-1 


dAp^ 


dAp^ _ 


-9 


dPa 


dPr P 


9-1 


dAp^ 


dAr 


1 


dT 


dPr P 


9-1 


dAr, 
ON 


dAM _ ^ 

dr] 





(3.19) 



The system above have one possible and convenient solution. 



1 

1" 

A„ = S = 0. 



At = t^Pt, Ap^ = 



Therefore, from the one-form canonical momenta above the new NC first-order Lagrangian 
can be computed directly. However, we also have to consider that the model remains 
second-order in velocities. Consequently, we will use the canonical momenta given by 

Arj = S . 

The respective NC first-order Lagrangian is. 



where the symplectic potential is written as 

V(0 = ^-3ka + Aa^ + PTa-^". (3.23) 

12a 

In order to obtain a commutative first-order Lagrangian, we propose a coordinate 

transformation in the classical phase space, analogous to the shift-operator Xj = Xi + ^9ijp' 
used in NC Quantum Mechanics (NCQM), given by 

Pa = Pt = t5f- (3-24) 
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Applying the transformation above in (3.22), we obtain the new commutative first-order 
Lagrangian, 

C = Pah + PtT + E?) - V{i) , (3.25) 
where the symplectic potential 

_p2 

y(0 = -j^ - ^ka + Ao^ + Pra"^" (3.26) 

is written in NC coordinates that satisfy the usual Poisson brackets given in Eq. (3.17). It is 
important to note that despite the variables in the Lagrangian Eq. (3.25) are commutative, 
there are NC contributions into the symplectic potential. Finally, the modified super- 
Hamiltonian of the system is identified as being the symplectic potential, so it can be 
written as _ 

H = -^^!^^^^^-3ka + Aa^ + PTa-^''. (3.27) 
12a 

Notice that when 6 = P = from (3.17) and (3.24), we recover, from (3.27), the Hamilto- 
nian in (3.3). 

4 Classical behavior of the NC cosmological models 

In order to investigate the contributions coming from the noncommutativity between the 
canonical variables and momenta in the classical FRW cosmological models, we derive the 
dynamical equations by computing the Hamilton's equations for the total Hamiltonian 
NH, where N is the lapse function and H is the modified super-Hamiltonian Eq. (3.27). 
We also use the constraint equation H = 0. The new momenta Pa and Pt, present 
in the H, are given by Eq. (3.24). In those expressions of Pa and Pt, and /3 are 
the parameters associated with the noncommutativity among the canonical variables and 
momenta, respectively (Eq. (3.17)). In the present study, we are going to consider only 
the contribution coming from the parameter f3. In other words, we shall fix 6 = 0. The 
reason is because we do not want the resulting dynamical equations depending on terms 
of order greater than two in the velocities. 

Assuming that = 1, we obtain the following two dynamical equations for the scalar 
factor a, 



a2 + a2 3+3 3 ' ^ ^ 

- + ^ + ^ = A-p-f3a--'-'d-a), (4.2) 
a 3 

where the dot means derivative with respect to the time coordinate t, in the present gauge. 

The first equation is the generalization for the NC models of the Priedmann equation. Both 

equations reduce the corresponding dynamical equations of the commutative models when 

we use that /3 = 0. 

It will be very useful to rewrite the generalized Friedmann equation (4.1) with the 
following form, 
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+ V{a) = 0, (4.3) 

where 

y(a) = A;--Aa2 + ^a-3"--a-3"-i, (4.4) 
3 3 3 

where C is a positive integration constant which is related to the initial fluid energy density 
(po). We notice that the total energy of this conservative system is equal to zero. From the 
observation of the potential curve V{a), we shall be able to derive the qualitative dynamical 
behavior of a{t). 

It is important to notice that the NC models described by the two dynamical equations 
above satisfy the energy conservation equation. In order to show this result, we use both 
Eqs. (4.1) and (4.2) to obtain, 

p + 3-{p + P) = 0. (4.5) 

a 

This is the energy conservation equation to the commutative version of the models. There- 
fore, the noncommutativity does not violate the energy conservation law. 

Both equations (4.1) and (4.2) are not independent. On the other hand, it can be 
shown that the set of solutions of Eq. (4.2) is not the same set of solutions of Eq. (4.1). 
However, the set of solutions of Eq. (4.2) is the same set of solutions of 



SnGp A_P _3a-2 
a2 ^ a2 3 3 3 

where f{t) = foa~^{t) and /o is an integration constant. We may impose that the solutions 
of Eq. (4.2) be, also, solutions of Eq. (4.1). It can be accomplished by fixing the following 
initial condition on the velocity a{t = to) = clq, 



««=^n-^+|+^-f««-'"-^K^ (4.7) 



where a{t = to) = oq. The only way this initial condition satisfy Eq. (4.6) is when 
/o = 0. Since this result must be valid for all times this initial condition guarantees that 
all solutions a{t) of Eq. (4.2) are also solutions of Eq. (4.1). In the following analysis we 
shall restrict our attention to the positive sign in Eq. (4.7). 

In order to derive the scalar factor as a function of t, we shall, initially, observe the 
potential curve V{a) from Eq. (4.1) and then solve Eq. (4.2). Unfortunately, for generic 
values of the different parameters present in Eq. (4.2), this equation does not have algebraic 
solutions. Therefore, we shall solve it numerically using a Fehlberg fourth-fifth order Runge- 
Kutta method with degree four interpolant, for each different values of the parameters. 

We have four parameters in Eq. (4.2). The first one is k, which is associated with 
the curvature of the spatial sections and we may assume three different values: -1, 0, +1. 
The second one is A, the cosmological constant, which can be positive, negative or zero. 
The third one is /3, the NC parameter, which can be also positive, negative or zero. It 
is important to mention that this last case /3 = means that the model is commutative. 
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Finally, we have the parameter a, which is present in the equation of state for the perfect 
fluid {p = ap). Each value of a defines a different perfect fluid. Here, we shall consider 
six different values of a = (1, 1/3, 0, —1/3, —2/3, —1), which represents respectively: stiff 
matter, radiation, dust, cosmic strings, domain walls and dark energy. Taking into account 
all possible values of the parameters, we have considered 162 qualitatively different cases 
and we solved Eq. (4.2) for all of them. In fact, we have solved Eq. (4.2) for a number 
of cases greater than 162 because we have considered not only the different signs for the 
cases of A and /3 but also for different absolute values of these parameters. 

After solving Eq. (4.2) for all possible cases, mentioned above, we have reached the 
following general conclusions. If the parameter (3 is positive, it has the net effect of pro- 
ducing a negative acceleration in the scalar factor evolution. In cases where the universe 
is expanding, the presence of a positive /?, will slow the expansion or even stops it and 
forces the universe to contract. If the parameter /3 is negative, it has the net effect of pro- 
ducing a positive acceleration in the scalar factor evolution. In cases where the universe is 
expanding, the presence of a negative /3, will increase the expansion speed. On the other 
hand, in cases where the Universe is contracting, the presence of a negative /3 may force the 
Universe to expand. From these observations, we notice that /3 modifies the acceleration 
of the scale factor in the opposite way that the cosmological constant does. Next, we shall 
present some particular cases where the above conclusions can be clearly verified. 

4.1 Universe with A: = 0, a = and A > 

For this case the Universe has flat spatial sections, it is filled with dust and the cosmological 
constant is positive. The potential V{a) in Eq. (4.4) is given by. 



From the potential expression it is easy to see that this universe starts to expand from 
a singularity at a = 0. For, 



the Universe expands initially in a decelerated rate and later on in an accelerated rate 
toward a oo. Since, in the present case C, A > 0, the above condition is satisfied for all 
cases where /? < 0. The condition Eq. (4.9) guarantees that the potential V{a) Eq. (4.8) is 
always negative. For /3 = 0, we have the commutative solution. If (3 is positive and satisfies 
Eq. (4.9), the universe expands in a smaller rate than in the commutative case. On the 
other hand, if /3 is negative the universe expands in a greater rate than in the commutative 
case. As an example we show in figure 1 the potential V{a), Eq. (4.8). For the case where 
A = 0.01, C = 0.01 and j3 = O.Ol(red), O(Wue), -Om{green). Then, in figure 2, we show 
the solutions to Eq. (4.2) for the present case with the same values of the parameters of 
figure 1. 





/5 < (27^2^/4) (1/=^) 



(4.9) 
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0,5^ _ I _ 1,5 _ 2 _ 2,5 _ 3 




Figure 1. V{a) Eq. (4.8), for C = 0.01, A = 0.01 and /3 = O.Ol(red), 0(6/ue), -0.01(5reen). 




Figure 2. Solutions a{t) to Eq. (4.2), for C = 0.01, fc = 0, A = 0.01, a = and 
(3 = Q.01{red),0{blue),~0.01{green). 
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4.2 Universe with k = 0, a = and A < 

For this case the universe has flat spatial sections. It is filled with dust and the cosmological 
constant is negative. The potential V{a) Eq. (4.4) is given by, 



where |A| is the cosmological constant absolute value. 

From the expression for the potential (4.10) it is easy to see that this universe starts to 
expand from a singularity at a = 0. Then, for all possible values of (3 the scale factor reaches 
a maximum value and starts to contract toward a big crunch singularity at a = 0. For /3 = 0, 
we have the commutative solution. If f3 is positive the scale factor is forced to contract in a 
stronger way than in the commutative case. Here, the maximum value of a is smaller than 
in the commutative case. On the other hand, if /3 is negative the scale factor will contract 
in a weaker way than in the commutative case. Here, the maximum value of a is greater 
than in the commutative case. As an example we show in figure 3 the potential V{a), Eq. 
(4.10), for the case where A = -0.01, C = 0.01 and /3 = O.Ol(red), O(Wne), -0.01{green). 
Then, in figure 4, we show the solutions to Eq. (4.2) for the present case with the same 
values of the parameters of figure 3. 




(4.10) 



Noncommutative cosmological models coupled. 



17 




Figure 3. V{a) Eq. (4.10), for C = 0.01, A = -0.01 and p = O.Ol(red), 0(6Zue), -O.Ol(green). 




Figure 4. Solutions a{t) to Eq. (4.2), for C = 0.01, k = 

O.Ol{red), 0{blue), ~0.01{green). 



0, A = -0.01, a = and /3 = 
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4.3 Universe with k = 0, a = 1/3 and A = 

For this case the Universe has flat spatial sections, it is filled with radiation and there is 
no cosmological constant. The potential V{a) Eq. (4.4) is given by, 

^(«) = -^ + f- (4-11) 
oa^ 6a 

Here, we have a very interesting situation. For /? < the universe starts to expand 
from a singularity at a = 0. It expands in a decelerated rate until it reaches asymptotically 
a —7- oo. For /? = 0, we have the commutative solution. This commutative case is well-know 
in the literature and it corresponds to Friedmann equation Eq.(4.3) that can be solved to 
give the algebraic solution: a(t) = At^^^'^\ where A is a constant. If /3 is negative the 
universe expands in a rate greater than the commutative case. On the other hand, if /? 
is positive the universe expands up to a maximum size given by c//3, then it is forced 
to contract toward a big crunch singularity at a = 0. Therefore, we see that the NC 
parameter may change drastically the universe evolution. As an example we show in 
figure 5 the potential V{a), Eq. (4.11), for the case where A = —0.01, C = 0.01 and 
P = 0. 01 {red),0 {blue), —0.01{green). Then, in figure 6, we show the solutions to Eq. (4.2) 
for the present case with the same values of the parameters of figure 5. It is important 
to mention that the term depending on (3 in Eq. (4.2) is not present in this case. The 
solutions of Eq. (4.2) depend on /3, in the present case, due to the initial condition Eq. 
(4.7). 
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-0,10" 



Figure 5. V{a) Eq. (4.11), for C = 0.01 and /3 = O.Ol(red), 0(6/Me), -0.01(5reen). 




Figure 6. Solutions a{t) to Eq. (4.2), for C = 0.01, fc = 0, A = 0, a = 1/3 and 
(3 = Q.Ol{red),0{blue), -O.Ol(green). 
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4.4 Universe with k = —1, a = —2/3 and A < 

For this case the Universe has constant negatively curved spatial sections, it is filled with 
a domain wall perfect fluid and there is a negative cosmological constant. The potential 
V{a) in Eq. (4.4) is given by, 

Via) = -1 - ^a' + ^a' - |a. (4.12) 

Here, we have the situation opposite to the previous case. For /3 > A, the universe 
starts to expand from a = 0, which is not a singularity. Then, the scale factor reaches a 
maximum value, 

amax = {C + VC^ + 12(/3 - A))/[2(/3 - A)] , 

and starts to contract toward a = 0. 

Since A < 0, the Universe is bounded for the commutative solution. If /3 is positive 
the scale factor is forced to contract in a stronger way than in the commutative case. Here, 
the maximum value of a is smaller than in the commutative case. 

On the other hand, if /3 < A the Universe starts to expand from a = and accelerates 
its expansion toward a — t- cx). Therefore, we will see another example where the NC 
parameter may change drastically the universe evolution. As an example we show in 
figure 7 the potential V{a), Eq. (4.12), for the case where A = —0.01, C = 0.01 and 
(3 = 0.011{red),0{blue), —0.011{green). Here, we can find an algebraic solution for the 
scalar factor as a function of time. Using Eqs. (4.1) and (4.2), and using the values of the 
parameters for the present case, we obtain the following differential equation, 

(J 

a + u'^a=—. (4.13) 
6 

The general solution of this equation is given by, 

a{t) = Ae''^' + Be-''^* + -^. (4.14) 

where A, B are integration constants and = (/3 — A)/3. 

If > Eq. (4.13) represents a driven harmonic oscillator with a constant driven 
force. The solution, Eq. (4.14), is oscillatory. On the other hand, if < Eq. (4.13) 
represents an unbounded system and the solution Eq. (4.14) grows exponentially. 
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Figure 7. V{a) Eq. (4.12), for A = -0.01, C = 0.01 and /3 = O.Ol(red), O(Wue), -~ 0.01 {green). 
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5 The noncommutativity and the Universe evolution 

Taking into account the above results we can see that the presence of /? in the cosmological 
equations may modify in a fundamental way the evolution of the Universe. In order to 
gain a better insight about the role of /3, let us rewrite the generalized Priedmann equation 
(4.1) in the following way, 



We may interpret A as a function of the scale factor that generalizes the cosmological 
constant. It is important to notice that, the contributions coming from A and the term 
depending on /3 oppose to each other. In order for them to reinforce each other, A and /? 
must have opposite signs. 

Depending on the value of a, we may have two entirely different scenarios for the 
evolution of the Universe described by A. 

5.1 First scenario 

For a > —2/3, the modulus of the term j3a~^'^~'^ decreases when a increases. This behavior 
agrees with the idea that noncommutativity should be most important at the beginning 
of the Universe. After some time, when the universe has evolved considerably a is large 
enough so that A ~ A/3. From this moment the Universe evolves without any remembrance 
of a NC phase. 

5.2 Second scenario 

For a < —2/3, the modulus of the term /3a~^"~^ increases when a increases. After 
some time, when the Universe has evolved considerably a is large enough so that A ~ 
— /3a~'^°~^/3. From this moment the Universe evolution is dominated by the NC term. 

If /3 > 0, the Universe will eventually reach a maximum size and contract again toward 
a — >■ 0. On the other hand, if /3 < 0, the Universe will expand forever in an accelerated 
rate. Here, we have a very interesting possibility. It is possible to consider the NC effects 
as a candidate to explain the present accelerated expansion of the Universe [43], [44]. In 
fact, some authors have already considered this possibility using classical, NC cosmological 
models different from ours [34], [35], [36]. 

As a matter of completeness, we observe that for a = —2/3, A is a constant. It is 
given by, A = (1/3) (A — /?). In this case, the NC parameter can be compared with the 
cosmological constant. If A = 0, the NC parameter plays the role of a cosmological constant 
with opposite sign. 




k SttGp 



+ A(a) 



(5.1) 



where 



A(a) = (l/3)(A-/3o- 



-3a-2 
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6 Final considerations and conclusions 

One of the mysteries that defies the theoretical physics today is how to promote the uni- 
fication of two independent sets of very different concepts. Namely, one set that rules the 
microscopic world, i.e., the quantum mechanics and the other that govern the macroscopic 
Universe, i.e., general relativity. The main motivation to unify both disconnected (so far) 
sets is to understand the physics of the early Universe and the physics that rules the ori- 
gin and the physical structure of the black holes. In few words, we are looking for the 
generalized idea of the so called quantum gravity. 

As we know the NC parameter has its value measured at the Planck scale. Besides, the 
string theory, a candidate to carry out this unification, embedded in a magnetic background 
(to explain it in a nutshell) , was found to have a NC algebra. In view of these two theoretical 
facts it is natural to believe that to investigate NC theories can be one of the adequate 
paths to conduct us to such unification. 

It is the main objective of this work to pursue this target. To accomplish this, we used 
the generalized NC symplectic formalism to introduce naturally noncommutativity inside 
the equations that provide the dynamics of the Universe, i.e., the Friedmann equations. 
With this procedure, we introduced a Planckian object inside the classical equations of 
motion of the Universe. 

We have constructed NC versions of several Friedmann- Robertson- Walker cosmological 
models. These models may have positive, negative or zero spatial sections curvatures. They 
are coupled to different types of perfect fiuids and they may have a positive, negative or 
zero cosmological constant. We have used the Schutz variational formalism in order to 
write the proper Hamiltonian for the models. We have performed a complete analysis of 
the evolution of the Universes described by each model. We solved the dynamical equations 
for the scale factor. 

Based on the results obtained here we have concluded that, if the parameter /? is 
positive, it has the final effect of producing a negative acceleration in the scalar factor 
evolution. In cases where the Universe is expanding, the presence of a positive /?, will slow 
the expansion or even stops it and forces the Universe to contract. If the parameter /? is 
negative, it has the effect of producing a positive acceleration in the scalar factor evolution. 
In cases where the Universe is expanding, the presence of a negative /3, will increase the 
expansion speed. On the other hand, in cases where the Universe is contracting, the 
presence of a negative /3 may force the Universe to expand. From these observations, 
we notice that /3 modifies the acceleration of the scale factor in the opposite way that 
the cosmological constant does. We have also presented some particular cases where the 
above conclusions could be clearly verified. Finally, based on our results, we construct a 
function (A) that depends on a, (3 and a and generalizes the cosmological constant. We 
have discussed the possible scenarios for the Universe evolution predicted by A, depending 
on the value of a. 
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